Nonlinear Dual-Mode Control of Variable- Speed 
Wind Turbines with Doubly Fed Induction 

Generators 

Choon Yik Tang, Member, IEEE, Yi Guo, Student Member, IEEE, and John N. Jiang, Senior Member, IEEE 



Abstract — This paper presents a feedback/feedforward non- 
linear controller for variable-speed wind turbines with doubly 
fed induction generators. By appropriately adjusting the rotor 
voltages and the blade pitch angle, the controller simultaneously 
enables: (a) control of the active power in both the maximum 
power tracking and power regulation modes, (b) seamless switch- 
ing between the two modes, and (c) control of the reactive power 
so that a desirable power factor is maintained. Unlike many 
existing designs, the controller is developed based on original, 
nonlinear, electromechanically-coupled models of wind turbines, 
without attempting approximate linearization. Its development 
consists of three steps: (i) employ feedback linearization to exactly 
cancel some of the nonlinearities and perform arbitrary pole 
placement, (ii) design a speed controller that makes the rotor 
angular velocity track a desired reference whenever possible, and 
(iii) introduce a Lyapunov-like function and present a gradient- 
based approach for minimizing this function. The effectiveness 
of the controller is demonstrated through simulation of a wind 
turbine operating under several scenarios. 

Index Terms — Wind energy, wind turbine, active power, reac- 
tive power, maximum power tracking, power regulation, nonlin- 
ear control. 

I. Introduction 

WIND power is gaining ever-increasing attention in 
recent years as a clean, safe, and renewable energy 
source. With the fast growth of wind generation in power 
systems, wind power is becoming a significant portion of the 
generation portfolio in the United States as well as many coun- 
tries in Europe and Asia 1 1 1 . Indeed, wind power penetration 
is planned to surpass 20% of the United States' total energy 
production by 2030 — a figure that is way beyond the current 
level of less than 5% |2|. Hence, to realize this vision, it is 
necessary to develop large-scale wind farms that effectively 
produce electric power from wind, and integrate them with 
the power systems. 

The integration of large-scale wind farms into a power 
system, however, changes the fundamental principle of its 
operation, which is to maintain reliability by balancing load 
variation with "controllable" generation resources. When a 
portion of these resources comes from "uncontrollable" wind 
generation, that portion of the resources can hardly be guar- 
anteed due to the intermittency of wind. As a result, the 
power system may fail to achieve the required balance. When 
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the level of wind power penetration is small, this issue may 
be safely neglected. However, with the anticipated increase 
in penetration, the issue becomes critical for power system 
reliability. As a case in point, half of the European grids 
experienced a severe difficulty in 2006 because several large 
wind farms, operating in the maximum power tracking (MPT) 
mode, produced excessive power that destabilized the grids. 
This event took place even though the wind penetration level, 
at that time, was low at only 7% |3|. 

Such an experience suggests that it is not advisable — 
and may be even disastrous — for a large-scale wind farm 
connected to a grid to always operate in the MPT mode, 
making wind turbines in the farm harvest as much wind energy 
as they possibly could, following the "let it be when the wind 
blows" philosophy of operation. Instead, it is highly desirable 
that the wind farm can also operate in the power regulation 
(PR) mode, whereby its total power output from the wind 
turbines is closely regulated at some desired setpoint, despite 
the fluctuating wind conditions. 

The ability to operate in the PR mode in addition to the 
MPT mode, as well as the ability to seamlessly switch between 
the two, offers many important advantages: not only does the 
PR mode provide a cushion to absorb the impact of wind 
fluctuations on total power output through power regulation, it 
also enables a power system to effectively respond to changes 
in reliability conditions and economic signals. For instance, 
when a sudden drop in load occurs, the power system may 
ask the wind farm to switch from the MPT to the PR mode 
and generate less power, rather than rely on expensive down- 
regulation generation. As another example, the PR mode, when 
properly designed, allows the power output of a wind farm to 
smoothly and accurately follow system dispatch requests, thus 
reducing its reliance on ancillary services such as reliability 
reserves. 

To enable large-scale wind farms to operate well in these 
two modes and switch seamlessly between them, numer- 
ous challenges must be overcome. This paper is devoted to 
addressing a subset of these challenges, by presenting an 
integrated framework for controlling the rotor voltages and 
the blade pitch angle of variable-speed wind turbines with 
doubly fed induction generators (DFIGs). The paper presents a 
feedback/feedforward nonlinear controller developed based on 
original, nonlinear, and electromechanically-coupled models 
of wind turbines, without attempting approximate lineariza- 
tion. The controller simultaneously enables: (a) control of the 
active power in both the MPT and PR modes, (b) seamless 
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switching between the two modes, and (c) control of the 
reactive power so that a desirable power factor is ensured. 
Its development consists of three steps. First, we show that, 
although dynamics of a wind turbine are highly nonlinear 
and electromechanically coupled, they offer a structure, which 
makes the electrical part feedback linearizable, so that arbitrary 
pole placement can be carried out. Second, we show that 
because the electrical dynamics can be made very fast, it is 
possible to perform model order reduction, so that only the 
first-order mechanical dynamics remain to be considered. For 
this reduced first-order model, a speed controller is designed, 
which enables the rotor angular velocity to track a desired ref- 
erence whenever possible. Finally, we introduce a Lyapunov- 
like function that measures the difference between the actual 
and desired powers and present a gradient-based approach for 
minimizing this function. The effectiveness of the controller is 
demonstrated through simulation of a wind turbine operating 
under a changing wind speed, changing desired power outputs, 
modeling errors, and noisy measurements. 

To date, a significant amount of research has been performed 
on the control of variable-speed wind turbines |4il- |i37J . The 
existing publications, however, are substantially different from 
our work in the following aspects: 

(i) The mechanical and electrical parts of the wind turbines 
are considered separately in most of the current litera- 
ture: P1- ll2n considered only the mechanical part, while 
Il22l - lf36l considered only the electrical part, focusing 
mostly on the DFIGs. In contrast, in this paper we con- 
sider both the mechanical and electrical parts. Although 
ll37l also considered both these parts, its controller 
was designed to maximize wind energy conversion, as 
opposed to achieving power regulation (i.e., only operate 
in the MPT mode). In comparison, our controller can 
operate in both the MPT and PR modes as well as 
seamlessly switch between the two. 

(ii) For those references Pi1- ll2n considering only the me- 
chanical part, control of the active power has been the 
main focus: ||4l, 131 maximized the wind power capture 
at low to medium wind speeds by adjusting both the 
generator torque and the blade pitch angle, fS^-f^l did 
the same by adjusting only the generator torque, and 
lfT3l - ll2l1 aimed to maintain the rated rotor speed and 
limit the power production at high wind speeds by 
controlling the blade pitch angle. Although the results 
are interesting, there is a lack of discussion on reactive 
power control, which is sometimes crucial to power 
system reliability. On the other hand, for those refer- 
ences Il22l - lf36l considering only the electrical part, $22\ . 
II23I considered only the active power control, while 
II24I - II27I dealt with the reactive power In contrast, 
we provide an integrated solution to the more difficult 
problem of simultaneously controlling both the active 
and reactive powers by appropriately adjusting both the 
rotor voltages and the blade pitch angle. 

(iii) From a modeling point of view, the wind turbine model 
we consider is one of the most comprehensive. Specifi- 
cally, [iT2l - lfT4l assumed a linearized model of the me- 



chanical torque from wind turbines, while H- lfTTI . ifTSi 
considered a nonlinear one. In addition, ll23l - lf32l . ll35l - 
||37J considered a reduced-order DFIG model, while 
II33I . Il34l considered a full-order one. In contrast, we 
consider a fifth-order, nonlinear, electromechanically- 
coupled model and attempt no linearization around some 
operating points, 
(iv) Finally, from a controls point of view, the control 
techniques we use are different from those adopted in 
B1- ll37l . For control of the mechanical part of the wind 
turbines, several techniques have been considered, in- 
cluding proportional-integral-derivative (PID) [12J-[14J, 
fuzzy control E), S-lHOl, El], adaptive control 
Q, robust nonlinear control 15], ifTS] , sliding mode 
control ifTTl . and wind-model-based predictive control 
ifTTl . Similarly, for control of the electrical part, various 
ways of controlling the DFIGs have been proposed, 
including vector/decoupling control with or without a 
position encoder (|23]-||25], E], 1321, IMl-lES], direct 
power control f28l, power error vector control f3U\, 
robust coordinated control |26|, linear quadratic regu- 
lator II29I . nonlinear inverse system method ll27l . and 
passivity -based control f33l|. In comparison, the con- 
troller we propose here consists of a mixture of several 
linear/nonlinear control techniques, including feedback 
linearization, pole placement, and gradient-based opti- 
mization and involving three time scales. The reason 
our controller is more complex is that the problem we 
address is inherently more challenging (see (i) and (ii)) 
and the model is more comprehensive (see (i) and (iii)). 
As it follows from the above brief review of the current 
literature, this paper contributes significantly to the state of 
the art on the control of wind turbines. 

The remainder of the paper is organized as follows: 
Section describes a model of variable-speed wind tur- 
bines with DFIGs. Section |lll] introduces the proposed feed- 
back/feedforward nonlinear controller. Simulation results are 
shown in Section |IV] Finally, Section |V] concludes the paper 

II. Modeling 

Consider a variable-speed wind turbine consisting of a 
doubly fed induction generator (DFIG) and a power electronics 
converter, as shown in Figure[T] The DFIG may be regarded as 
a slip-ring induction machine, whose stator winding is directly 
connected to the grid, and whose rotor winding is connected 
to the grid through a bidirectional frequency converter using 
back-to-back PWM voltage-source converters. 

The dynamics of the electrical part of the wind turbine are 
represented by a fourth-order state space model, constructed 
using the synchronously rotating reference frame (dg-frame), 
where the relation between the three phase quantities and the 
dq components is defined by Park's transformation 1.381 . The 
voltage equations are [39l 

Vds = Rsids - l^sVqs + '^'■Pdsi (1) 
d 

Vqs = Rslqs + ^^s^ds + -TlVqs, (2) 
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Fig. 1. Schematic of a typical variable-speed wind turbine with a DFIG. 



Vdr = Rridr - (t^s " i^r)'^qr + '^^dr, (3) 

, . d 

Vqr = Rrlqr + (t^s " ^rWdr + '^^<ir^ W 

where Vds, Vqs, Vdr, Vqr are the d- and g-axis of the stator and 
rotor voltages; ids, iqs, idr, iqr are the d- and g-axis of the 
stator and rotor currents; ipds, ^qs, fdr, fqr are the d- and 
g-axis of the stator and rotor fluxes; ujs is the constant angular 
velocity of the synchronously rotating reference frame; lu^ is 
the rotor angular velocity; and R^, Rr are the stator and rotor 
resistances. The flux equations are ||39l 

(5) 
(6) 
(7) 
(8) 



fds = Lgids + L 

m'^dr: 

^qs ^s^qs ^m^qr, 
^dr ^m^ds ^~ ^r^dr^ 



where Lg, L^, and L,„ are the stator, rotor, and mutual 
inductances, respectively, satisfying Lg > L,n and Lj. > L^- 
From (l5]l-(|8]l, the current equations can be written as 

Ids = —rVds — T-rVdr, (9) 

(JLjs (TLgLir 



'~Pqri 



1 



-^ds 



-^dr, 



Lm. 1 

-fqs ^ T^qr, 



(7 L a Lr 



(10) 

(11) 

(12) 



where cr = (1 — ^ '£ ) is the leak coefficient. Selecting the 
fluxes as state variables and substituting (l9]l-(fT2]i into ([rii-(|4|i, 
the electrical dynamics in state space form can be written as 



d 


Rs 

l-'fids 

aLg 


+ ^s^qs 


RsLm 
CtL s 


d 


-LOgipds - 


Rg 


1 Rs L'fn 
CtL g Lr 


d 


RfLfn 
J J Vds 


Rr , , 
-;—Vdr + (Ws 


d 

dt"^"^^ 


(7 L g Lf 


- (Ws - 


^r)Vdr — 



Vqr 



Rr 
(jLr 



^qs, 



(14) 



(15) 



uqr ■ 



(16) 



Neglecting power losses associated with the stator and rotor 
resistances, the active and reactive stator and rotor powers are 



given by f^Ol 



Qs — '^qs'^ds '^ds'^qs-) 
Pr ^ ^'^dr'^dr ~ '^qr'^qri 
Qr " '^qr'^dr ~t~ '^dr'^qr-i 



(17) 
(18) 
(19) 
(20) 



and the total active and reactive powers of the turbine are 

P^Pg+Pr, 



(21) 
(22) 



where positive (negative) values of P and Q mean that the 
turbine injects power into (draws power from) the grid. 

The dynamics of the mechanical part of the wind turbine 
are represented by a first-order model 



1± -T 

dt 



fUJr, 



(23) 



where C/ is the friction coefficient, T,„ is the mechanical 
torque generated, and is the electromagnetic torque given 
by 



^qs'^ds ^ds'^qs i 



(24) 



where positive (negative) values mean the turbine acts as a 
generator (motor). The mechanical power captured by the wind 
turbine is given by pT| 



P — T 



3 

w ' 



(25) 



where p is the air density; A = ttR^ is the area swept by the 
rotor blades of radius R; Vw is the wind speed; and Cp{X, jS) is 
the performance coefficient of the wind turbine, whose value 
is a function |41 1 of the tip speed ratio A, defined as 



A 



(26) 



as well as the blade pitch angle j3, assumed to lie within some 
mechanical limits Pmin and /3max- This function is typically 
provided by turbine manufacturers and may vary greatly from 
one turbine to another 141]. Therefore, to make the results of 
this paper broadly applicable to a wide variety of turbines, 
no specific expression of Cp{X, (3) will be assumed, until it is 
absolutely necessary in Section IIVI to carry out simulations. 
Instead, Cp(A, /3) will only be assumed to satisfy the following 
mild conditions: 

(Al) Function Cp(A, f3) is continuously differentiable in both 
A and /3 over A € (0,oo) and /3 e [/3min, /3max]- 

(A2) There exists c € (0, oo) such that for all A G (0, oo) 
and 13 e [Pmin, Pmax], we have Cp(A, /3) < cA. This 
condition is mild because it is equivalent to saying that 
the mechanical torque T„i is bounded from above, since 
r,„ oc according to ^ and (|26|. 

(A3) For each fixed (3 e [/5min, /3max], there exists Ai E 
(0, oo) such that for aU A € (0, Ai), we have Cp(A, 13) > 
0. This condition is also mild because turbines are 
designed to capture wind power over a wide range of 
A, including times when A is small. 
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For the purpose of simulation in Section |IV| the following 
model of Cp(A, /3) will be assumed, which is presented in 
and also adopted in MATLAB/SimuUnk R2007a: 



and (fT3Tl-(fT6]l in matrix forms as follows: 



Cp(A, /3) = ci(^^ - C3/3 - 04^6 +ceX, 



where 



1 

A" 



1 



0.035 



A + 0.08/3 /33 + r 



(27) 



(28) 



T 
T 



and the coefficients are ci — 0.5176, C2 = 116, C3 — 0.4, 
C4 = 5, C5 = 21, and cg = 0.0068. 

As it follows from the above, the wind turbine studied in this 
paper is described by a fifth-order nonlinear dynamical system 

with states [(pds fqs fdr fqr ^r]'^ , COntrols [Vdr Vqr /3]"^, 

outputs [P Q]^, "disturbance" Vw, nonlinear state equations 
(fT3Tl-(fT6b and ( |23] ), and nonlinear output equations (fT7])-(l22b. 
Notice that the system dynamics are strongly coupled: the 
"mechanical" state variable w,- affects the electrical dynam- 
ics bilinearly via ( fTSI ) and ( fTSI l. while the "electrical" state 
variables [tpds fqs fdr ^qrV affect the mechanical dynamics 
quadratically via (l9b-(fT2b. (|24| ) and ( l23b . Since the stator wind- 
ing of the DFIG is directly connected to the grid, for reliability 
reasons [vds Vqs]^ are assumed to be fixed, i.e., not to be 
controlled, in the rest of this paper Moreover, since (l9]l- (fT2] | 
represent a bijective mapping between Ypds 'Pqs fdr fqr 
and [ids iqs idr iqrV and since the currents [ids iqs idr iqr] 
the rotor angular velocity lOj., and the wind speed can all be 
measured, a controller for this system has access to its entire 
states (i.e., full state feedback is available) and its disturbance 
(i.e., the wind speed Vw)- A block diagram of this system is 
shown on the right-hand side of Figure |2l 



III. Controller Design 

In this section, a feedback/feedforward nonlinear controller 
of the form depicted on the left-hand side of Figure |2] is 
presented. By adjusting the rotor voltages Vdr and u^,. and the 
blade pitch angle /?, the controller attempts to make the active 
and reactive powers P and Q track, as closely as possible — 
limited only by wind strength — some desired, time-varying 
references Pd and Qd, presumably provided by a wind farm 
operator. When Pd is set to sufficiently large, i.e., larger than 
what the turbine can possibly convert from wind, it means 
the operator wants the turbine to operate in the MPT mode; 
otherwise, the PR mode is sought. The value of Qd, along with 
that of Pd, reflects a desired power factor PF^ = , ^ the 

operator wants the turbine to also maintain. 

The controller development consists of three steps, which 
are described in Sections IIII-Al - |lircl respectively. 



A. Feedback Linearization and Pole Placement 

For convenience, let us introduce the variable x 
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(29) 



X, 



(30) 



where A is a constant matrix and, as was pointed out at the 
end of Section [III both Vds and Vqs are constants not to be 
controlled. Note that the only nonhnearities in ( |29] l are the 
two products of the state variables, i.e., —01^X4 and w^xa. 
Also note that these nonhnearities appear on the same rows as 
the control variables Vdr and Vqr- Thus, feedback linearization 
(431 may be used to cancel them and subsequently perform 
arbitrary pole placement (HI, i.e., let 



Vdr 



'Jqr 



UJrXi — 
-UJrXs 



Kfx 



Ul, 
Ko X + U2, 



(31) 
(32) 



[Xi Xi X3 



YPdt 



'fqs fdr "ySgr]"^ and rewrite (l9Tl-(fT2]i 



where Ki,K2 G M , the first terms on the right-hand side 
of ( I3T] ) and ( l32b are intended to cancel the nonhnearities, the 
second terms are for pole placement, and the third are new 
control variables ui and U2, to be designed later. 

To implement dSTT i and ( |32] |. full state feedback on the fluxes 
X and the rotor angular velocity lOj. are needed. While the latter 
is relatively easy to measure, the former is not. Fortunately, 
this difficulty can be circumvented by first measuring the 
currents — which is feasible — and then calculating the fluxes 
from (Is])-®. This explains the fourth input of the nonlinear 
controller block in Figure |2l 

Substituting (EB and (O into (|29l) yields 

X ^ {A - BK)x + [vds Vqs Ul U2]'^ , (33) 

where B [02x2 -^2x2]^ and K = [Ki K2Y is the state 
feedback gain matrix. Since the electrical elements in the 
DFIG are physically allowed to have much faster responses 
than their mechanical counterparts, K in ( |33] | may be chosen 
so that A — BK is asymptotically stable with very fast eigen- 
values. With this choice of K and with relatively slow-varying 
Ul and U2 (recall that Vds and Vqs are constants), the fourth- 
order hnear differential equation ( l33T l may be approximated 
by the following static, linear equation: 

x = -{A- BK)-^[vdsVqsUiU2]^ . (34) 

As a result, the fifth-order model described in (fT3]l-(fT6]l and 
( I23I ) may be approximated by the first-order model described 
in ( I23] ) along with ( l34l ). As will be shown next, this approxi- 
mation greatly simplifies the design of ui and U2- Therefore, 
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Fig. 2. Structure of the multivariable, feedback/feedforward nonlinear controller, developed based on original, nonlinear dynamics of the wind turbine. 



we will assume, in the sequel, that K is chosen so that the 
electrical dynamics ( |33]) are asymptotically stable and so fast 
that they may be approximated by ( l34b . 

B. Tracking of Desired Angular Velocity 

The second step of the controller development involves con- 
structing a speed controller that ensures the angular velocity 
of the rotor, ujr, tracks a desired, time-varying reference, LOrd, 
whenever possible. The construction may be divided into four 
substeps as described below. 

Substep 1. First, we show that the electromagnetic torque 
Te defined in may be expressed as a quadratic function 
of the new control variables ui and U2. From (l24l) and 
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(35) 



Equation ( |35] | suggests that Tg is a quadratic function of x, 
while (l34l i suggests that x, in turn, is an affine function of ui 
and U2, since Vds and Vqs in (l34l i are assumed to be constants. 
Hence, Tg must be a quadratic function of ui and U2- Indeed, 
an explicit expression can be obtained as follows: since A — 
BK is asymptotically stable and thus nonsingular, it may be 
written as 
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where each dij depends on A, B, and K. From (|34lt-(|36]>. 
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where qi, q2, 93, bi, 62, and a are constants defined as 
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92 = o _r r (.di3di4 + ^14^43 - 0^23^34 - ^24^33), (39) 

(40) 
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Substep 2. Next, we show that the quadratic function ( l37b 
relating ui and U2 to Tg has a desirable feature: its associated 
Hessian matrix [ql ql] is always positive definite, regardless of 
the parameters of the electrical part of the DFIG, as well as 
the choice of the state feedback gain matrix K. The following 
lemma formally states and proves this assertion: 

Lemma 1. The Hessian matrix [^293] ( 1371 ) is positive 
definite. 

Proof: From (O, and gOl), the 

determinant of A — BK and the leading principal minors 91 
and 51(73 — q2 of [ql ql] can be written as 



1^ - BK\ 



91 

91 93 - 92 



d2 r4 
2 _ ^s^m 



(44) 
(45) 
(46) 



where 



A = ( - RsL„iki2 + {LsLr - Ll^)ki3 - RsLrku 

+ RrL,+R,Lr) Ai 
+ {RsLmkll + RsLrki3 + {LsLr - im)fcl4 

+ RsRr-LsLr + L^^)A2, (47) 
Ai = RsLmk2l + RsLrk23 + {LsLr - im)fc24 
+ RrLs + RsLr, 
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A2 — RsLmk22 + { — LgLr 
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and kij is the zj entry of K. Since A — BK is nonsingular, 
is > L„i, and > im, it follows from ( l44l i that A 7^ 0. 
Since A 7^ 0, it follows from (|46] l that qiqz — > and 
from ( |47| i that Ai and A2 cannot be zero simultaneously. The 
latter, along with ( l45T l, implies that gi > 0. Since qi > and 
(Zi^a ~ 92 > 0, [ga 93 ] in < l37] i must be positive definite. ■ 
Substep 3. Next, we show that there is a redundancy 
in the control variables ui and U2, which may be exposed 
via a coordinate change. Observe from (l23T l that the first- 
order dynamics of ujr are driven by Tg. Also observe from 
Substeps 1 and 2 that is a convex quadratic function of 
ui and U2- Thus, we have fwo coupled control inputs (i.e., 
ui and U2) affecting one state variable (i.e., a;^), implying 
that there is a redundancy in the control inputs, which may 
be exploited elsewhere (to be discussed in Section IIII-Cl l. To 
expose this redundancy, first notice that because the Hessian 
is positive definite, it can be diagonalized, i.e.. 



matrix 



92 93 



there exist an orthogonal matrix M containing its eigenvectors 
and a diagonal matrix D containing its eigenvalues, such that 



qi 92 

92 93 



M = D. 



(48) 



Indeed, 

M = 
where 



A2-g3 



22 



+ ^/{>^2-q3Y+ql, 



D 



Ai 





A2 



Al,2 — 



91 + 93 ± V(9i +93)^ -4(9193 -ql) 



Next, consider the following coordinate change, which trans- 
forms Ml S M and U2 € M in a Cartesian coordinate system 
into r > and 6' G [0, 2tt) in a polar coordinate system: 



zi' 



where 



1 



(49) 



(50) 



In terms of the new coordinates r and 9, it follows from (|37] | 
and (|48ll-(|50ll that 



where 



a' = a - - [61 62] 



(51) 



9i 


92 


-1 




92 


93. 







Using ([38Tl-(|43]l, a' may be simplified to 

^UJsRs 

implying that it is always negative. Comparing dSTT l with dJTj i 
shows that the coordinate change ( |49] l and ( ISOl l allows us to de- 
couple the control variables, so that in the new coordinates, r is 
responsible for driving the first-order dynamics of ujr through 



Te of ( fSTT l, while 6 does not at all affect (and, hence, is 
redundant as far as the dynamics of uj^ are concerned). The 
design of r and 9 will be discussed in Substep 4 and Section 
IIII-CI respectively. 

Substep 4. Finally, a speed controller is presented, which 
ensures that the rotor angular velocity tUr tracks a desired 
time-varying reference uJrd, to be determined in Section ITlI-CI 
provided that Urd is not exceedingly large. Combining 
and ( BTl i yields 



C 



where, according to (IZST l. 



fUJr, 



Tm{uJr,l3,V^) 



(53) 



(54) 



Here, Tm is written as Tm{(^r, f3,Vw) to emphasize its de- 
pendence on LUr, P, and Vw Observe from (|53] | that, if the 
control input were reaZ-vaZMec/ instead of being nonnegative, 
feedback hnearization may be applied to cancel all the terms 
on the right-hand side of (l53T l and insert linear dynamics 
a{LiJr — <-^rd)^ i-S-. we may let 



= T„i{uJr,l3,Vw) - a' - CfUJr + a{uJr - w^d), 

so that 



(55) 



(56) 



By letting the controller parameter a be positive, ( |56] | implies 
that ujr always attempts to go to uJrd- Unfortunately, because 
cannot be negative, the speed controller ( fSST l — and, hence. 



the linear dynamics (1561 1 — cannot be realized whenever the 
right-hand side of ( ISSl l is negative. To alleviate this issue, (l55l l 
is slightly modified by setting to zero whenever that occurs, 
i.e., 

r ^ = max{rm (w^ , ^5, K; ) - a' - C'/w^ + a (w^ - i:^rd) , 0} . 

(57) 

Notice that (ISTl l contains a feedforward action involving the 
"disturbance", i.e., the wind speed Vw This explains the first 
input of the nonlinear controller block in Figure |2] 

To analyze the behavior of the speed controller dSTT i. sup- 
pose UJrd, P, and Vuj are constants and consider the function 
g, defined as 

g{uJr, l3, Vnj) = TrniuJr, /?, K,) - fl' - CfUJr- (58) 

The following lemma says that 9(0;^, /3, Ku), when viewed 
as a function of ujr, has a positive root ujr^\ below which 
g{uJr, f3,Vn]) is positive: 

Lemma 2. For each fixed (3 € [/3min, /3max] '^nd T4j > 0, 
there exists ivi^^ G (0,cxd) such that g(iLii^\ /S^Vw) = and 
g{ujr, 13, Vw) > for all ujr € (0, uji^''). 



Proof: Due to the fact that a' in (|52] | is negative, there 
exists ojr.i such that —a' — CfUr > for all Ur E (0,Wr,i). 
Due to Assumption (A3) of Section lUl (l26]l, and (Elll, there 
exists ijjr.2 such that T,n{iOr, (3, Ku) > for all ujr G (0,^^,2)- 
Hence, from (|58] |. we have 5(0;^, /?, Ku) > for all £ 
(0, minjoj^.i, '^r.2})- In addition, due to Assumption (A2), 
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(|26] l, (|54] i. (|58T l. and a' < 0, there exists ujr.3, sufficiently 
large, such that g{ujr.3, f3,Vw) < 0. These two properties of 
g, along with Assumption (Al) and the Intermediate Value 
Theorem, imply that there exists at least one positive root ujr 
satisfying g{ujr: P.V^) = 0. Letting Wr^-* be the first of such 
roots completes the proof. ■ 
The following theorem, derived based on Lemma |2] says 
that as long as the desired rotor angular velocity ojrd is not 
exceedingly large, i.e., does not exceed the first root Lui^^ of 
g{iL!r, /?, Vw), the closed-loop dynamics ( |53] | and (|57] | have an 
asymptotically stable equilibrium point at ujrd- 

Theorem 1. Consider the first-order dynamics (I53l l and the 

speed controller (157b . Suppose LUrd, cind Vu, are constants, 
with LOrd satisfying < Uj-d < w^^-*. Then, for all Wr(O) > 0, 

Proof: Substituting ( |57] | into (l53T l and using dSST l yield 



Critical root u-v'^' 



Jiur = min{a(wrd 

,(1) 



(59) 



Suppose < uJrd < LOr . We first show that ujr = i^rd is the 
unique equilibrium point of (|59] |. Suppose = w^d. Then, 
a{ujrd — oJr) in (|59] | is zero, whereas ^(wr,/?, Ku) in (l59] l 
is positive, due to Lemma |2] Thus, cj^ ~ 0, implying that 
is an equilibrium point. Next, suppose < ujr < ^rd- 
Then, a{ujrd — ^r) is positive, and so is /?, Ku), due 

again to Lemma |2] Hence, ujr > 0, implying that there 
is no equilibrium point to the left of Urd- Finally, suppose 
bJr > ^rd- Then, a{ujrd — i^r) is negative. Therefore, Ur < 0, 
implying that there is no equilibrium point to the right of 
UJrd- From the above analysis, we see that tUr = ujrd is the 
unique equilibrium point of (|59] l. Next, we show that the 
equilibrium point ujr — ujrd is asymptotically stable in that 
for all ujr{0) > 0, limf^oo u!r{t) — ujrd- Consider a quadratic 
Lyapunov function candidate V : (0,oo) — > R, defined as 

1 



y{^r) = -{uJr- i^rd) , 



(60) 



which is positive definite with respect to the shifted origin 
UJr = UJrd- From ( |59] l and ( |60] l. 



V{uJr) 



UJrd) T[im{a{uJrd - UJr),giuJr, (3, V^)) . 

(61) 



Note that whenever < < '-^rd^ ct{^rd ~ '^r) > and 
g{uJr, f3,Vw) > 0, so that V{ujr) < according to (|6T| i. On 
the other hand, whenever ujr > ujrd, ct{^rd ~ i^r) < 0, so that 
V{ujr) < 0. Finally, when ujr — ujrd, V{ujr) = 0. Therefore, 
V{ujr) is negative definite with respect to the shifted origin 
UJr = UJrd- It follows from [43 1 that ujr = uird is asymptotically 
stable, i.e., for all ujr{0) > 0, limf^oo w^ft) = ujrd- ■ 
Theorem [T] says that the first root a;^^ is a critical root, 
for which utrd should never exceed, if we want ujr{t) to go 
to UJrd regardless of ujr{0)- Figure |3] shows, for the MAT- 
LAB/SimuHnk R2007a model of Cp(A,^) given in ^ and 
(l28T l. how the critical root wf^' depends on (3 and V^- Notice 
from the figure that Wr^"* is insensitive to /3 but proportional 
to Vw, meaning that the larger the wind speed, the higher the 
"ceiling" on the desired rotor angular velocity. Also notice that 




3540 1 ' ' ' ' ' ' ' ' ' ' ' 1 7080 

2 4 6 8 10 12 14 16 18 20 22 24 
Wind speed Vw {m/s} 

Fig. 3. Critical root lj'^' as a function of blade pitch angle /3 and wind 
speed Vw 



Wr^-* of more than 3500 in the per-unit system is extremely 
large, meaning that for this particular turbine there is no need 
to be concerned about w^d exceeding uji^\ 

C. Lyapunov-like Function and Gradient-based Approach 

The third and final step of the controller development in- 
volves introducing a Lyapunov-like function, which measures 
the difference between the actual and desired powers, and 
utilizing a gradient-based approach, which minimizes this 
function. 

Recall from the beginning of Section that the objective 
of the controller is to make the active and reactive powers, P 
and Q, track some desired references, Pd and Qd, as closely 
as possible. In the MPT mode, where the goal is to generate as 
much active power as possible while maintaining an acceptable 
power factor, Pd is set to a value that far exceeds what the 
wind turbine can possibly produce (e.g., in the per-unit system, 
Pd > 1), while Qd is set to a value representing the desired 
power factor PF^ — . . In this mode, making P and 

Q approach Pd and Qd is equivalent to maximizing the active 
power output while preserving the power factor. In the PR 
mode, where the goal is to regulate the powers, both Pd and 
Qd are set to values representing power demands from the 
grid. In this mode, making P and Q approach Pd and Qd 
amounts to achieving power regulation. Hence, the values of 
Pd and Qd reflect the mode the wind farm operator wants the 
wind turbine to operate in. However, as far as the controller 
is concerned, it does not distinguish between the two modes; 
all it does is try its best to drive P and Q to Pd and Qd- 

To mathematically describe the aforementioned controller 
objective, consider the following positive definite, quadratic 
Lyapunov-like function V of the differences P — Pd and Q — 
Qd: 



V = l[P-Pd 



Wp Wpq 




P - Pd 


Wpq Wq 




Q - Qd_ 



(62) 



>o 



where Wp, Wq, and Wpq are design parameters that allow one 
to specify how the differences P — Pd and Q — Qd, as well as 
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their correlation {P — Pd){Q ~ Qd), should be penalized. With 
this V, the above controller objective can be restated simply 
as: make V go to zero, because when this happens, P and Q 
must both go to Pd and Qd- Since it is not always possible 
to achieve this objective — due to the fact that the wind may 
not always be strong enough — below we will attempt instead 
to make V as small as possible by minimizing it. 

To minimize V, we first show that F is a function of uJrd, 
6, P, Vu,, Pd, and Qd, i.e., 

V ^ fiu;rd,0,l3,V^,,Pd,Qd) (63) 

for some /. Note from ( |62] | that V depends on P, Q, Pd, and 
Qd- Also note from (fT7]i-(l20li. (IHT i. and ^ that P and Q, 
in turn, depend on i, Vdr, and Vqr (recall that Vds and Vqs are 
constants). Thus, 

V = fi{i,Vdr,Vqr,Pd,Qd) (64) 

for some fi. Next, note from (l9Tl-(fT2Ti that i depends on x; 
from (ISTT i and ( |32] | that Vdr and Vqr depend on x, ujr, wi, and 
U2', and from (|34] | that x further depends on ui and U2- Hence, 

{i,Vdr,Vqr) = h{^r,Ui,U2) (65) 

for some /2. Furthermore, note from (|49] | and (fSOl l that wi and 
it2 depend on r and 9, where r, in turn, depends on oj^, w^d, 
/3, and Vw through dSTT i. Therefore, 



(66) 



for some /a. Finally, assuming that w^d does not exceed the 
first root Wr^'' and assuming that ujrd, P, and T4, are all rel- 
atively slow-varying (see below for a discussion). Theorem [T] 
says that ujr goes to uJrd- Thus, after a short transient. 



UJr 



^rd- 



(67) 



Combining (I64ll-(l67]i. ( l63T l is obtained as claimed. 

Now observe that the first three variables [ujrd, S, (3) in 
(l63T l are yet to be determined, while the last three variables 
{Vw, Pd,Qd) are exogenous but known. Therefore, for each 
given {Vw, Pd,Qd), (wrd,6',/3) can be chosen correspond- 
ingly in order to minimize V. This defines a mapping from 
{Vni,Pd,Qd) to {ujrd,0,P), i.e., 

iurd,0,(3)^F{V,,,Pd,Qd) 

= argmin^^^^^^^^^) /(xi,X2,X3,K,,Fd,Qd)- (68) 



In principle, the mapping F in (I68t may be constructed 
either analytically, by setting the gradient of /(•) to zero and 
solving for the minimizer {ujrd, S, (3) in terms of {Vw, Pd, Qd), 
or numerically, by means of a three-dimensional lookup table. 
Unfortunately, the former is difficult to carry out, since /, 
being composed of several nonlinear transformations (|64]|- 
(I&TI i. has a rather complex expression. On the other hand, the 
latter is costly to generate and can easily become obsolete due 
to variations in system parameters. More important, selecting 
{uJrd,d,f3) as a static function of {Vw,Pd,Qd) as in (l68T l 
may lead to steep jumps in {LUrd,0,P) because V^ is ever- 
changing and may change dramatically, and both Pd and Qd 
from the wind farm operator may experience step changes. 
Such steep jumps are undesirable because large fluctuations 
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i^rd 
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'^drj ^qr 



Fig. 4. Internal structure of the proposed nonlinear controller. 

in UJrd may prevent Ur from tracking it, while discontinuous 
changes in (3 may be mechanically impossible to realize, cause 
intolerable vibrations, and substantially cut short the lifetime 
of the turbine blades. 

To alleviate the aforementioned deficiencies of selecting 
{u!rd,6,f3) according to ( |68] ), a gradient-based approach is 
considered for updating (w^d, 6*, /3): 

df 



^rd 



(2 



-£3 



du]rd 

dl 

9/3' 



(69) 
(70) 
(71) 



where £1,62, £3 > are design parameters, which are meant 
to be relatively small, especially ei and £3, in order to avoid 
steep changes in Urd and /3. The partial derivatives g^p^, fg, 
and II in dlDi-dZB can be calculated in a straightforward 
manner using (|64]|-(|67]|. but are omitted from this paper due 
to space limitations. These partial derivatives are practically 
implementable since, like /, they depend on Lo^d, S, (3, V^, 
Pd, and Qd, all of which are known. With this gradient-based 
approach, {uJrd, 6, (3) is guaranteed to asymptotically converge 
to a local minimum when {Vw,Pd,Qd) is constant, and track 
a local minimum when {Vw,Pd,Qd) varies. 

To help the readers better understand the proposed nonlinear 
controller depicted in Figure |2] and described in this section, 
the internal structure of this controller is revealed in Figure |4] 
Observe that each arrow in this figure represents a signal, 
whereas each tiny box represents equations relating the signals. 

IV. Simulation Studies 

To demonstrate the effectiveness of the controller presented 
above, MATLAB simulations have been carried out, in which 
the controller is applied to a 1.5 MW GE turbine with 575 V 
base voltage and 60 Hz base frequency. To describe settings 
and results of the simulations, both the per-unit system and the 
physical unit system will be used, given that they are popular 
in the literature. 

The simulation settings are as follows: All values of 
the wind turbine parameters are taken from the Wind Tur- 
bine Block of the Distributed Resources Library in MAT- 
LAB/Simulink R2007a. Specifically, the values are: a;s(pu) = 
1, i?s(pu) = 0.00706, i?r(pu) = 0.005, L,(pu) = 3.071, 
i^(pu) = 3.056, Lm(pu) = 2.9, J(pu) = 10.08, C/(pu) = 
0.01, Anon, = 8.1, and Cp_nom = 0.48. In addition, the 
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Fig. 5. Scenario 1 illustrating the maximum power tracking (MPT) mode. 



base wind speed is Vwj,-^^^ — 12m/s. In MATLAB/Simulink 
R2007a, the mechanical power captured is given by 

p„(pu) ^ ^w.nd_basePno:. ^^^^^^ ^^^^^^3^ (72) 
-'elec_base 

where Pwmd_base = 1.5 MW, Pnom = 0.73, and Pelec_base = 

1.5 X lO'^/O.gVA. From it can be seen that at the 

base wind speed Vu,(pu) = 1, Pm(pu) is capped at 0.657. 
Furthermore, without loss of generality, the constant stator 
voltages are set to i>ds(pu) ~ 1 and ^^^(pu) = 0. 

For the proposed controller, we let the desired poles of the 
electrical dynamics (|33] | be located at —15, —5, and —10 ± 
5j, so that the coiTesponding state feedback gain matrix K, 
calculated using MATLAB's place () function, is 

_ r 5135.9 259.2 20.3 1.9" 
~ [-2676.7 4289.9 -1.3 19.7 ' 

In addition, we let Wp = 10, Wq = 1, and Wpq = 0, implying 
that we penalize the difference between P and P^j much more 
than we do Q and Q^. Finally, we choose the rest of the 
controller parameters as follows: a = 10, ei = 4 x 10^'^, 
£2 1 X IQ-*, and eg = 2. 



Based on the above wind turbine and controller parameters, 
simulations have been carried out for four different scenarios. 
Description of each scenario, along with the simulation result, 
is given below: 

Scenario 1: Maximum power tracking (MPT) mode. In this 
scenario, we simulate the situation where the wind speed 
Vw experiences step changes between 12m/s and 7.2 m/s, 
while the desired powers Pd and Qd are kept constant at 
1.5 MW and 0.15 MW, so that the desired power factor is 
PFd = 0.995. Since P„ cannot exceed 0.657 x 1.5 MW at 
the base wind speed Vw_bass = 12 m/s, the wind turbine is 
expected to operate in the MPT mode. Figure |5] shows the 
simulation result for this scenario, where the key signals are 
plotted as functions of time in both the per-unit and physical 
unit systems wherever applicable. Observe that, after a short 
transient, the wind turbine converts as much wind energy 
to electric energy as it possibly could, as indicated by Cp 
approaching its maximum value of 0.48 in subplot 2 (which 
translates into P approaching its maximum possible value in 
subplot 3). Also observe that, when Ku goes from 12 m/s to 
7.2 m/s and from 7.2 m/s back to 12 m/s, Cp drops sharply 
but quickly returns to its maximum value. Note from subplot 
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Fig. 6. Scenario 2 illustrating the power regulation (PR) mode. 



4 that, regardless of Vw, the power factor PF is maintained 
near the desired level of 0.995. Moreover, note from subplots 

5 and 6 that the angular velocity tracks the desired time- 
varying reference uj^d closely (subplot 6 is a zoom-in version 
of subplot 5). Finally, the control inputs Vdr, Vqr, and (3 are 
shown in subplots 7 and 8, respectively. Note that, to maximize 
Cp, (3 is kept at its minimum value jiaun — deg. 

Scenario 2: Power regulation (PR) mode. In this scenario, 
we simulate the situation where is kept constant at the 
base value of 12m/s, while Pd experiences step changes from 
0.45 MW to 0.3 MW and then to 0.6 MW, and Qd is such that 
PFd = 0.995. Since Pd is always less than 0.657 x 1.5 MW at 
the base wind speed of 12 m/s, the wind turbine is expected to 
operate in the PR mode with different setpoints Pd- Figure |6] 
shows the simulation result for this scenario. Observe from 
subplot 2 that Cp is less than its maximum value of 0.48. 
This suggests that the wind turbine attempts to capture less 
power than what it possibly could from wind, since Pd is 
relatively small. Indeed, as can be seen from subplots 3 and 
4, the turbine produces just enough active and reactive powers, 
making P track Pd closely while maintaining PF at PF^. Also 
observe from subplots 5 and 6 that lu^ closely follows uj^d, as 
desired. Finally, note from subplot 8 that /3 increases slightly 



in order to capture less power between 1200s and 2400s, when 
Pd is smallest. 

Scenario 3: Seamless switching between the MPT and PR 
modes. In this scenario, we simulate the situation where Pd ex- 
periences large step changes between 1.5 MW and 0.75 MW, 
Qd again is such that PF^ is 0.995, and an actual wind profile 
from a wind farm located in northwest Oklahoma is used to 
define V^. The actual wind profile consists of 145 samples, 
taken at the rate of one sample per 10 minutes, over a 24- 
hour period. In order to use this wind profile in a 1-hour 
simulation (as in Scenarios 1 and 2), we compress the time 
scale, assuming that the samples were taken over a 1-hour 
period. Note that compressing the time scale in this way makes 
the problem more challenging because the wind speed varies 
faster than it actually does. Figure |7] shows the simulation 
result for this scenario, with subplot 1 displaying the wind 
profile. Observe from subplots 2 and 3 that, for the first 1200 
seconds during which Pd is 1.5 MW, the turbine operates in 
the MPT mode, grabbing as much wind energy as it possibly 
could, by driving Cp to 0.48 and maximizing P. At time 
1200s when Pd abruptly drops from 1.5 MW to 0.75 MW, the 
turbine seamlessly switches from the MPT mode to the PR 
mode, quickly reducing Cp, accurately regulating P around 
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Fig. 7. Scenario 3 illustrating the seamless switching between the MPT and PR modes under an actual wind profile from a wind farm located in northwest 
Oklahoma. 



Pd, and effectively rejecting the "disturbance" Vw Note that 
between 2100s and 2400s, the wind is not strong enough to 
sustain the PR mode. As a resuh, the MPT mode resumes 
seamlessly, as indicated by Cp returning immediately to its 
maximum value of 0.48. Finally, at time 2400s when Pd goes 
from 0.75 MW back to 1.5 MW, the turbine keeps working in 
the MPT mode, continuing to maximize both Cp and P. Notice 
from subplots 4-6 that, over the course of the simulation, both 
PF and are maintained at PFd and w^d, respectively, despite 
the random wind fluctuations. Also notice from subplot 8 that 
/3 increases somewhat during the PR mode in order to help 
capture less power. 

Scenario 4: Robustness of the proposed controller. In this 
scenario, we simulate the exact same situation as that of 
Scenario 3 (i.e., with the same Vw, Pd, and Qd) but with 
modeling errors and measurement noise. That is, we allow 
for modeling errors in the friction coefficient C/ and the 
performance coefficient Cp (due, for example, to changing 
weather conditions, blade erosions, and aging) as well as 
measurement noise in the wind speed (since is usually 
measured by an anemometer located on the naceUe behind 
the blades of a wind turbine). Specifically, we assume that 
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Fig. 8. Contour plots of the nominal and actual Cp(A, /9) for Scenario 4. 



the nominal Cf used by the controller is 0.01 (pu), whereas 
the actual C/ used in the simulation is 0.012(pu), so that 
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Fig. 9. Scenario 4 illustrating the robustness of the proposed controller to modeling errors in C/ and Cp and noisy measurements in Vw 



Cf has a 20% modeling error. Moreover, we assume that the 
nominal Cp used by the controller is given by ( |27| i and ( |28] l 
with ci = 0.5176, C2 = 116, = 0.4, C4 = 5, C5 21, and 
ce = 0.0068, whereas the actual Cp used in the simulation is 
also given by (|27]i and (|28]l but with ci = 0.45, C2 = 115, 
C3 = 0.5, C4 = 4.5, C5 = 22, and cg = 0.003. Figure [8] 
displays the contour plots of the nominal and actual Cp{X, jS) 
for A e [2, 15] and /3 e [0, 15], showing that Cp has noticeable 
modeling errors. In particular, the nominal Cp attains its 
maximum of 0.48 at (A,/3) = (8.1,0), whereas the actual Cp 
attains its maximum of 0.39 at (A, /3) — (8.45, 0). Finally, we 
assume that the measured Vw used by the controller, denoted 
as Vu,_meas, IS related to the actual Vw used in the simulation 
via 

K;_meas(i) = K; W + 0.5 + 0.5 sin(0.5i) + 0.25cos(t), 

where the second term on the right-hand side represents 
a constant measurement bias, while the third and fourth 
represent measurement noises with different amplitudes and 
frequencies. Figure |9] shows the simulation result for this 
scenario. Comparing this figure with Figure |7] the following 
observations can be made: first, Cp in Figure |9] attains its 
maximum value of 0.39 in the MPT mode, as opposed to the 



0.48 attained by Cp in Figure Q Second, PF in Figure |9] has a 
larger fluctuation compared to PF in Figure |7] but nonetheless 
is maintained around PF^;. Third, oj^ in Figure|9]does not track 
ojj-d as closely as uj^ in Figure|7]does. Nevertheless, despite the 
wind fluctuations, modeling errors, and noisy measurements, 
the controller performs reasonably well, as evident by how 
close Cp is to its maximum value of 0.39 in the MPT mode, 
how close P is to Pd in the PR mode, and how close PF is 
to PFd throughout the simulation. Therefore, the controller is 
fairly robust. 



As it follows from Figures |5]-|9] and the above discussions, 
the proposed controller exhibits excellent performance. Specif- 
ically, the controller works well in both the MPT mode under 
step changes in the wind speed (Scenario 1) and the PR mode 
under step changes in the power commands (Scenario 2). In 
addition, it is capable of seamlessly switching between the 
two modes in the presence of changing power commands 
and a realistic, fluctuating wind profile (Scenario 3). Finally, 
the controller is robust to small modeling eiTors and noisy 
measurements commonly encountered in practice (Scenario 4). 
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V. Conclusion 

In this paper, we have developed a feedback/feedforward 
nonUnear controller, which accounts for the nonlinearities of 
variable- speed wind turbines with doubly fed induction gen- 
erators, and bypasses the need for approximate linearization. 
Its development is based on applying a mixture of linear 
and nonlinear control design techniques on three time scales, 
including feedback linearization, pole placement, and gradient- 
based minimization of a Lyapunov-like potential function. 
Simulation results have shown that the proposed scheme not 
only effectively controls the active and reactive powers in both 
the MPT and PR modes, it also ensures seamless switching 
between the two modes. Therefore, the proposed controller 
may be recommended as a candidate for future wind turbine 
control. 
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